The bilinear-biquadratic spin-1 chain undergoing quadratic Zeeman effect 
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The Heisenberg model for spin-1 bosons in one dimension presents many different quantum phases 
including the famous topological Haldane phase. Here we study the robustness of such phases in 
front of a SU(2) symmetry breaking field as well as the emergence of unique phases. Previous 
studies have analyzed the effect of such uniaxial anisotropy in some restricted relevant points of the 
phase diagram. Here we extend those studies and present the complete phase diagram of the spin-1 
chain with uniaxial anisotropy. To this aim, we employ the density matrix renormalization group 
(DMRG) together with analytical approaches. The complete phase diagram can be realized using 
ultracold spinor gases in the Mott insulator regime under a quadratic Zeeman effect. 
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I. INTRODUCTION 

Quite generally, ultracold atoms trapped in optical lat- 
tices provide a toolbox to reproduce various kinds of 
Hubbard models 1 - offering thus a unique possibility to 
simulate quantum magnetism 2 ^—. As it is well known, 
Hubbard models reduce to various spin models in certain 
limits. Most of these limits, and even many others, are 
accessible with ultracold lattice atoms. 

There are several ways of achieving such a situation, 
the simplest perhaps is the case where spin- 1/2 is as- 
sociated to the presence or the absence of a structure- 
less boson in a lattice site (the, so called, hard-core 
boson limit). More generally, if one considers systems 
of particles with several internal states in many body 
Mott insulator states, such systems are well described by 
low energy effective Hamiltonians that incorporate super- 
exchange interactions^, and quite generally correspond to 
spin models 2 ^— , or composite fermions models^^. If the 
internal states correspond to pseudo-spin, the resulting 
Hamiltonians typically do not have special symmetries. 
If, on the contrary, the internal states correspond to the 
degenerate manifold of the hypcrfine atomic level F, i.e. 
{\F, tuf)} of an ultracold spinor gas, the resulting Hamil- 
tonians, in the absence of symmetry breaking fields, are 
SU(2) symmetric. The corresponding spin Hamiltoni- 
ans are given by powers of the nearest neighbor Heisen- 
berg interactions H = J2k a k(Si • Sj) fc (for F = 1 see 
for instance^ 1 -, for F = 2,3/2,5/2 c.frU, and for general 
Fermi systems c.fji^r— ). 

The are two kinds of symmetry breaking fields that 
can be realized and controlled both in atomic-ionic and 
in condensed matter systems: those corresponding to the 
linear Zeeman effect for atoms (z aligned magnetic field 
in condensed matter), or the quadratic Zeeman effect 
(single ion anisotropy in condensed matter). The linear 
Zeeman effect is important, but not so interesting, since 
the magnetization is a constant of the motion, and it 
can effectively be "gauged out"—. Also, in experiments, 



frequently magnetic fields oscillate rapidly in time, and 
their effect averages to zero. This is not the case for 
the quadratic Zeeman effect, which leads to more com- 
plex and interesting effects* 1 ^— Particulary interesting 
are these effects in one dimension (ID), where the role of 
quantum fluctuations is enhanced. 

For these reasons, considerable efforts have been de- 
voted to the study of the simplest case of spin-1 chains, 
namely the, so called, bilinear-biquadratic spin-1 chains 
undergoing a quadratic Zeeman effect. On the other 
hand, for the most of the known spin-1 alkali atomic sys- 
tems, the scattering lengths in different scattering chan- 
nels are practically equal, so that the scattering is ef- 
fectively spin state independent. That means that the 
systems exhibits higher symmetry, in the case of F = 1 
the SU(3) symmetry. A similar situation occurs with 
earth alkali fermionic atoms, that can have spin F and 
exhibit spin state independent collisions, i.e. SU(2F+ 1) 
symmetry^. 

Very recently, Rodriguez et alJ£ have studied the 
phase diagram of the bilinear-biquadratic spin-1 chains 
undergoing the quadratic Zeeman effect close to the 
SU(3) symmetric point. These authors have used an 
elegant effective field model, that allowed them to eval- 
uate approximately the phase boundaries for the case of 
vanishing magnetization in any dimension. The study of 
Refi 1 ^ was, however, restricted to a quite limited region 
of parameters. The reason is that they consider only the 
region of the phase space achievable in alkaline ultracold 
gases by means of standard Fcshbach resonances. 

The present paper complements and extends the study 
of Rodriguez et alJ£ to the full range of parameters of the 
uniaxial spin-1 chain. We point out that by using the idea 
of employing "quantum simulators in excited states"—, 
i.e. studying properties of the spin chain not necessarily 
in the ground state, but in the state of maximal energy, it 
is possible to cover the entire range of parameters for the 
bilinear biquadratic spin-1 chains: from ferromagnetic to 
fully anti-ferromagnetic. The idea is here very simple: 
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The state of maximal energy, if gapped, is dynamically 
stable, and can only undergo thermodynamic instability, 
which, however, requires interactions with a reservoir etc. 
The maximal energy state will thus remain stable if the 
system is well isolated. 

The phase diagram in the spin chain is obtained 
numerically in the entire parameter region using a fi- 
nite size density matrix renormalization group (DMRG) 
algorithm^. Using degenerate perturbation theory, we 
show in a neat way how the emergence of the phases can 
be obtained from the anisotropy of an effective XXZ spin- 
1/2 model. Our results are summarized in Fig. 1, where 
for reference we mark also the region studied ini£. Our 
results agree quantitatively with the description provided 
in — in this region, although we found some significant 
discrepancies in the boundaries of some of the phases, as 
we shall comment later. 

The paper is organized as follows: in Sec. [II] we in- 
troduce the model and the quantities of interest. Here 
we briefly review the earlier, mainly condensed matter, 
literature for a non-zero anisotropy field. In Sec. Ill we 
initiate the study of the complete phase diagram. We 
discuss the results in various important limits, first for 
large negative D (Sec. MI A[) where perturbation theory 
leads to an effective XXZ spin 1/2 model. The distinct 
emerging phases are subsequently analyzed: dimer phase 
(Sec. [TTTT3]) : the Neel phase (Sec. InTCj) ; the Haldane 
phase (Sec. IIII D[) ; the ferromagnetic phases (Sec. IIIIEI) 
and the critical phases (Sec. IIII F|) . We conclude in Sec. 

m 



II. THE SPIN-1 CHAIN IN A UNIAXIAL 
FIELD: PREVIOUS RESULTS 
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FIG. 1: (Color online) Top: Phase diagram of a spin-1 chain 
in a Zeeman quadratic field as a function of 8 and D obtained 
by means of the DMRG algorithm. The dashed red box indi- 
cates the portion of phase diagram studied in Ref.— . Bottom: 
the anisotropy A of the effective spin- 1/2 XXZ model in the 
limit of large negative I? as a function of 9 (see the results 
of perturbation theory in Sec. IIII A|) . The vertical dotted 
lines separate the resulting phases: Neel (A > 1), critical 
( — 1 < A < 1) and ferromagnetic (A < —1). 



We consider a one dimensional spin-1 chain of length 
L with open boundary conditions whose Hamiltonian is: 

H = H BB + Hu = 

L-l 

= JJ2 [cos(9)Si • Sj+i + sin(0)(S 4 • S t+1 ) 2 ] 
i=i 

L 

+ Dj^^ (1) 

i=l 

where Si = (S X i, S y i, S Z i) are the ith site angular mo- 
mentum operators. The model is the sum of two terms, 
Hqb and Hjj. The former, proportional to J, is known 
as the bilinear-biquadratic spin-1 Hamiltonian, while the 
latter is a uniaxial field term proportional to D. In the 
rest of this section we review known results about this 
model. 



A. The case D = 

For D = 0, the phase diagram of Hbb as a function 
of 8 G [— 7r; 7r] is well known (see for example^). For 



— 7T < 8 < — 37r/4 and for 7r/2 < 8 < tt the ground 
state is ferromagnetic with a net spontaneous magnetiza- 
tion along a direction which breaks spontaneously space 
isotropy. For — 37r/4 < 8 < —ir/4 the ground state is 
characterized by a non vanishing dimer order parameter: 

V=\H l -H t+1 \ (2) 

where H l = cos(<9)Sj • S l+1 + sm(8)(S l ■ S i+1 ) 2 . In this 
phase, translational invariance is broken and singlets of 
neighboring spins appear. On the other hand, a nematic 
phase has been predicted in the region between the dimer 
and ferromagnetic phased. Its existence is however still 
debated^—. The difficulty in establishing the true na- 
ture of the region close to the ferromagnetic phase is due 
to the fact that the dimer order parameter decays expo- 
nentially as 8 — > — 37r/4. 

For — it /A < 8 < 7r/4 the system is in the Haldane 
phase. This is a topological phase, containing for 8 = 
arctan(l/3) the Affleck, Lieb, Kennedy, Tasaki (AKLT) 
points whose ground state can be written in terms of 
matrix product states. The entire phase is character- 
ized by free edge spins and a non vanishing string order 
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III. THE COMPLETE PHASE DIAGRAM 



i+r — 1 



0= lim (S Z i exp[z7r > *S Z7 -]5 



(3) 



j=i+i 



Recently, it has been shown that topological phases in 
one dimensional systems can be detected by a doubly de- 
generate entanglement spectrum^. In other words, the 
eigenvalues of the reduced density matrix of a portion of 
the chain are always doubly degenerate in a topological 
phase. Moreover, the characterization of the phase tran- 
sitions around the Haldane phase in terms of the Schmidt 
gap, i.e. the difference of the two largest density matrix 
eigenvalues, has been recently reported in^i. 

Finally for tt/A < 9 < tt/2 the system is in a criti- 
cal, gapless phase, characterized by zero energy modes 
at momentum q = 0,±27r/3 and has been extensively 
studied 2 ^ The existence of these gapless modes is re- 
flected in a peak at q — 2tt/3 in the magnetic and nematic 
structure factors: 



Sc(q) 



cxp[iq(m - n)](C m C„ 



(4) 



where C = S z or C = S 12 respectively. The periodicity 
induced by these momenta is also connected to the ex- 
istence of trimers in the phase. Very recently, an order 
parameter has been proposed which can be experimen- 
tally observed in an optical lattice setting^ 2 -. 



B. 



The case D ^ 



The Hamiltonian of Eq. (|T|) has also been studied for 
9 = and different values of D22r— . The resulting 
ground state phase diagram contains three phases: in the 
vicinity of D = the system is in the already mentioned 
Haldane phase; for large negative values of D, the ground 
state exhibits a spontaneous staggered magnetization: 



(5) 



and is therefore in an Ising-Neel phase; finally, for large 
positive D the ground state is adiabatically connected to 
the trivial state 1 00 . . . 0) in which all the spins are in the 
S Z n = cigenstatc. For this reason this phase has been 
dubbed the "large D" phase. 

Recently, Rodriguez et alJ£ have proposed to real- 
ize model dJ) for D ^ with spinor ultracold atoms 
in optical lattices. They found that by making use of 
Feshbach resonances in 87 Rb atoms, the interval 9 E 
(—tt + arctan |; — 7r/2) around the SU(3) point can be 
realized. They studied the corresponding phase dia- 
gram in ID, 2D and 3D, finding unique phases: an XY- 
ferromagnetic phase for D > and 9 < —3n/4 and an 
XY-nematic phase for D < and 9 e [— 3-7T/4; — 7r/2]. 



In this paper we study the complete phase diagram of 
Hamiltonian (TT|) in the interval 9 G [— tt; tt] and D/J e 
[—2; 3]. Our main result is summarized in Fig. [1] showing 
the resulting phase diagram. In the limit of large nega- 
tive D/J we apply perturbation theory to reconstruct the 
phase diagram. In the general case, we conduct exten- 
sive numerical simulations using the finite size density 
matrix renormalization group (DMRG) algorithm (see 
for example 2 ^ and references therein) with open bound- 
ary conditions. The number of block states retained has 
been always chosen so that the discarded weight, con- 
nected to the algorithm accuracy, is less than 10~ 6 or 
until convergence is reached. Typically we used between 
100 and 250 states. In the remainder of the paper we dis- 
cuss the different phases individually and the transitions 
separating them. 



A. Limit of large negative D 

Let us start by considering the limit of large negative 
D so that we can apply perturbation theory assuming 
\J/D\ -C 1. For J = and negative D the ground state 
of Hjj is degenerate and pertains to the manifold spanned 
by the states: 



£1^2 ■■■e L ) 



(6) 



where e; = ±1. The unperturbed ground state energy 
is simply DL < where L is the size of the chain. In 
this limit, Hbb acts as a perturbation and we can apply 
first order degenerate perturbation theory. This consists 
in diagonalizing the perturbation Hamiltonian in the de- 
generate subspace. Projecting the Hamiltonian Hbb on 
this manifold leads to the effective spin- 1/2 Hamiltonian: 



H B B = -JLsin9 + 

i 

-(2cos0- sin 0)<t? ( 



isin#«a^ 1 +afa^ 1 )+ 



E (9)+H xxz (7) 



where <rf,crf and af are the spin- 1/2 Pauli matrices for 
the ith spin in the basis {|1) , | — 1)}. 

Therefore, apart from a constant energy shift Eq(6), for 
large negative D the system is described by an effective 
XX Z spin-1/2 model with anisotropy A = (2cos# — 
sin#)/| sin#|. As shown in Fig. [lj the quantity A(0) as 
a function of 9 spans twice the real axis when 9 changes 
from — 7r to tt. The properties of the XXZ models are 
well known^: it is ferromagnetic for A < —1, it is in 
the critical XY phase for — 1 < A < 1 and in the Ising- 
Neel phase for A > 1. We can, therefore, reconstruct 
the phase diagram in this regime as a function of 9 as 
sketched in Fig. [T] The system is in the Neel phase for 
— tt/2 < 9 < 7r/4, in the ferromagnetic phase for — n < 
9 < -37r/4 and for 7r/4 < 9 < tt and in the critical XY 
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FIG. 2: (Color online) (a) Finite size extrapolation of the 
dimer order paramter for 9 = — 0.4-7T and two values of D: 
D = — 0.2J (plus) for which the extrapolated result is Doo = 
0.27 and D = — 0.7J (crosses) for which the extrapolated 
value is negative and therefore X>oo = 0. Also shown are the 
corresponding fitting functions, (b) The extrapolated dimer 
order parameter for 6 = — 0.6-7T and large negative D. The 
plot is shown as a function of —D/J and in a log-log plot. 
The data are compared with a power-law with exponent 0.75 
(solid line). 



phase for -37r/4 < 9 < -n/2 and for 7r/4 < 9 < n/2. We 
use these analytical results as a guide for the numerical 
simulations in the limit of large negative D. Notice that 
in the opposite limit of large positive D, the ground state 
is non degenerate and given by the product state in which 
all the spins are in the state |0). In this case, perturbation 
theory gives poor results— In the following sections, 
we analyze separately each phase characterizing it with 
appropriate orders. 



With this technique we reconstruct the dimer phase as 
shown in Fig. [I] Surprisingly enough, the dimer phase 
persists even at D = —2 J. Although the dimer order pa- 
rameter is different from zero, it is really small as shown 
in Fig. [2jb) and decays not faster than algebraically to 
zero as D — > — oo. This is consistent with our predictions 
from perturbation theory that dimcrization is zero in this 
limit. We cannot however exclude that the border of the 
dimer phase occurs for a finite value of D/J < —2. This 
result does not agree with RefJ£ which for 9 = — O.bir 
predicts the transition around D = —0.3 J. 

Wc analyze the robustness of the dimer phase for posi- 
tive D (see the inset of Fig. [1]): for example for 9 = — 0.5n 
the dimer order parameter is already zero for D > 0.05 J. 
Again, this is in clear contrast with the results of RefJ£ 
which predicts the transition for D ~ 0.6J and a cusp- 
like behavior which we do not observe. The disagreement 
could be due to the small system sizes taken in Ref 
The finite size effects are in fact evident in Fig. Ufa), 
showing that a careful extrapolation to the thermody- 
namic limit is in order. 

Notice that the procedure we use here is reliable only 
away from criticality and serves us to approximately lo- 
cate the phase boundaries. To get more accurate results 
for the critical points as well as critical exponents we will 
make use in Sec. MI CI of finite size scaling theory. 

We leave as an open question the behavior of the 
dimer phase boundaries close to the point D = and 
9 = — 3/4-7T where the gap and the dimer order param- 
eter are expected to decrease exponentially. Accurate 
boundaries are therefore hard to get with our code and 
probably more sophisticated tensor networks techniques 
arc required. 



B. The dimer phase 

Wc begin our numerical analysis with the dimer phase. 
As said in Sec. Hi] this phase is characterized by the pres- 
ence of the dimer order parameter T>. In order to obtain 
results in the thermodynamic limit we computed T>, using 
DMRG, for different system sizes L = 96, 132, 168, 204 
and then performed a finite size extrapolation with a 
quadratic function: 



V 



L 



L 2 



(8) 



Two examples of such extrapolation arc shown in 
Fig. EJa) for 9 = -0.4tt: (i) for D = -0.2J the ex- 
trapolated dimer order is D M = 0.27J and therefore we 
conclude that this point is in the dimer phase; (ii) for 
D = —0.7 J the extrapolated result is zero and therefore 
the point is not in the dimer phase. We have compared 
our extrapolated result in the point D = and 9 = — it/ 2 
with the analytical result valid in the thermodynamic 
limit V = 1.1243 found in Ref4i (see also^). Our result 
is T>oo = 1.13 and the less than 1% relative error gives 
an idea of the quality of our numerical simulations. 



C. The Neel phase 

As found in previous works^— and predicted by 
our perturbation theory results, in the asymptotic limit 
D —> — oo, a Neel phase occurs. This phase is character- 
ized by a non zero spontaneous staggered magnetization 
defined in Eq. ([5]) which wc use as the order parameter. 
Similar to the dimer phase, to locate the boundaries of 
the Neel phase, we extrapolate the staggered magnetiza- 
tion to the thermodynamic limit. In order to compare our 
results to precise quantum Monte Carlo simulations^ we 
make use of finite size scaling (FSS) theory^ to obtain 
the corresponding critical exponents. In the so called 
scaling regime, the dependence of an observable Q, when 
a parameter g of the Hamiltonian (in our case 9 or D) 
varies close to the critical point g c , is given by: 



Q~L-W''f Q (\g-g c \L 1 / v 



(9) 



where v characterizes the divergence of the correlation 
length when approaching the critical point: £ ~ \g — 
g c \~ u and /3q is the order parameter exponent: Q ~ \g — 
g c \^ Q in the ordered phase. To locate the critical point 
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FIG. 3: (Color online) Finite size scaling of the staggered 
magnetization M z for lengths: L = 96, 132, 168, 204 for 9 = 
close to the critical point separating the Neel and Haldane 
phases. Inset: the corresponding collapse plot (see text). 



g c we plot QL ai for different sizes L as a function of g 
and change ol\ until there is a crossing of all the curves. 
This procedure will give us g c at the crossing point and 
Pq/v = ai. In order to compute /3q, we plot QL^ Q ^ U as 
a function of (g — g c )L a2 and change ol^ until we observe 
collapse of the data for all lengths L. In this way we 
extract j3 = a\jai and v = l/a^. 

Care must be taken, because in a finite size system the 
Zi symmetry of the Hamiltonian is not spontaneously 
broken. To circumvent this technical issue, we add a 
small magnetic field acting only on the first spin, corre- 
sponding to an additional term O.OlJS^i. 

We use this method for = so that we can com- 
pare our results with Refi^, which employs quantum 
Monte Carlo simulations and in which the estimate for 
the critical point is at D^ M = —0.316 J and for the expo- 
nent for the staggered magnetization pj^ 1 = 0.147 while 
v® M = 1.01 which are not far from the universal expo- 
nents of the 2D classical Ising model (/3a/, = 1/8 and 
v = 1). Our results for the staggered magnetization are 
shown in Fig. [3l Crossing and collapse of the data for 
lengths L = 96, 132, 168,204 happen for D c = -0.316 J, 
$ Mt = 0.11 ±0.01, v = 1.1 ±0.1 which agree with Ref^ 
and with the Ising universality class. We checked that for 
other fixed values of 8, the critical points, as D changes, 
are also of the Ising type. 

D. The Haldane phase 

The Haldane phase is the archetypal example of a topo- 
logical phase in one dimensional spin systems. In this 
paper we characterize it using two key quantities: the 
string order parameter defined in Eq. and the degen- 
eracy of the entanglement spectrum defined as follows. 



After finding the ground state \iPg) of the system using 
DMRG, we construct the reduced density matrix of half 
chain: 

Pl=Ti-r\iPg)(iPg\ (10) 

The entanglement spectrum is given by the set of eigen- 
values Ai of the density matrix pl such that ^ i Ai = 1 
and ordered in decreasing order. A topological phase, 
such as Haldane, is signaled by a doubly degenerate spec- 
trum. A handful quantity for detecting the Haldane 
phase is the Schmidt gap, given by the difference of the 
two largest eigenvalues of pl: AA = Ai — A2. 

We studied the robustness of the topological phase in 
presence of the SU(2) Zccman quadratic field using both 
the Schmidt gap and the string order parameter. In 
Fig. [4] we display our results for 9 = and for posi- 
tive and negative values of D/J. Both the Schmidt gap 
and the string order O clearly identify the Haldane phase, 
however one has to be careful since O is different from 
zero also in the Neel phase. This is evident if one con- 
siders the perfect Neel state |±1, — 1, ±1, — 1 . . .). This 
ambiguity is resolved by looking at the staggered magne- 
tization Mf 1 also reported in the plot. Contrary to the 
Neel phase, the string order goes instead to zero in the 
large D phase, and at the same time the Schmidt gap 
reopens as expected. Fig. [1] shows that the Schmidt gap 
and the staggered magnetization scale very similarly. In 
fact in a recent work^I we have shown that, in proxim- 
ity of the Neel-Haldane boundary, the two quantities Mf ' 
and AA scale with the same critical exponents and can be 
therefore considered to be two possible order parameters 
for the description of this transition. 

Finally, we have compared our numerical results ap- 
proaching the point 6 = 7r/4 and D = known as 
the Uimin-Lai-Sutherland (ULS) model^ with exact an- 
alytical calculations. The ULS model in the presence 
of a linear and a quadratic field was studied in^ and 
the two critical values of D are D + = 2.8284J and 
D- = — 0.9800J. Notice that these values have been 
obtained by dividing the corresponding results in^i by a 
factor a/2 because of the different Hamiltonian definition 
in our model. Our estimates arc D+ = 2.83 ± 0.05 and 
£)_ = —0.95 ±0.05 which are in excellent agreement with 
the analytical predictions. 



E. The Ferromagnetic phases 

The system is in a trivial ferromagnetic phase for D < 
and for — it < 6 < — Sir/A and 7r/2 < 6 < tt. In this Ising 
ferromagnetic phase, the model is characterized by a non 
zero magnetization M z = 1/ L^2 n S zn along z (except 
along the line D = 0, where the magnetization could be 
in any direction). 

On the contrary, for the same values of 9 but for 
D > 0, the net magnetization is zero along any direc- 
tion. However for small values of D the system is in 
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FIG. 4: (Color online) The Schmidt gap, the string order pa- 
rameter and the staggered magnetization for 9 = are shown 
as a function of D/J. The data have been obtained from a 
quadratic extrapolation for the lengths L = 96, 132, 168, 204 
(see Eq. l[8|). The shaded area represents the scaling regions 
close to the critical points where the quadratic extrapolation 
does not work. 



an AF-Ferromagnetic state as predicted in Ref.— . Al- 
though the "in plane" magnetization is zero, since the 
system is one dimensional, the correlations in the XY 
plane decay algebraically as shown in Fig. [5] To find the 
boundary between the XY-Ferromagnetic phase and the 
large D phase we study the magnetic structure factor for 
the XX correlations: 



S x (q) = 2^2cxp{iq(m 

mn 



^0] {^xrn$xn, 



(ii) 



In the XF-Ferromagnetic phase, the correlations are ex- 
pected to decay as (S xn S x (n+ r )) ~ r~ a . It is very hard 
to extract an accurate value of a by fitting the correla- 
tion functions. Instead, by considering that the magnetic 
structure factor diverges for q = with the system size: 
S x (q = 0) ~ one can better estimate a. In the 

large D phase, correlations decay exponentially with a 
finite correlation length and the corresponding magnetic 
structure factor S x (q = 0) is independent of L. There- 
fore, in order to distinguish the XY- Ferromagnetic phase 
from the large D phase, it is sufficient to study the be- 
havior of S x (q = 0) with L and check whether it diverges 
or not. An example of such analysis is shown in Fig. [5] 
for 6 = — 0.87T and D/J = 0.1. In the top panel the 
correlations {S xn S x ( n+r }) are shown for n = L/A as a 
function of r. The value of n has been chosen so as to 
minimize finite size effects from the boundaries. In the 
bottom panel, we show the magnetic structure factor at 
zero momentum as a function of the system size. Using 
a power-law fit we find a = 0.07. 
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FIG. 5: (Color online) Top: In-plane correlations 
{S X nS x ( n + r \) as a function of the spin distance r for L = 
168,0 = -0.87T, D/J = 0.1. Notice that r is limited to the 
interval between 1 and L/2 to avoid boundary effects. Inset 
the same plot in logarithmic scale. The data are compared 
with a powerlaw ~ 1/r 07 . Bottom: the magnetic structure 
factor for the XX correlations S x (q = 0) for zero momentum 
as a function of the system size L. The solid line as been 
obtained with a power-law fit ~ I/ 1- "' and the best-fit value 
is a = 0.07. 



F. The critical phases 

The phase diagram contains also two interesting criti- 
cal phases as shown in Fig. [1] denoted as Critical A and 
Critical B. We find interestingly enough that Critical A 
has at least two different central charges c ~ 1 and c ~ 2, 
each identifying a different universality class^. We base 
our results on the scaling of the von Neumann entropy 
with the size of a partition of the chain. If pi is the re- 
duced density matrix of a block of spins of length I when 
the global state of L spins is the ground state, then the 
von Neumann entropy is given by: 



S(l) = -Tcpi logp. 



(12) 



For critical points in one dimensional quantum spin 
chains, the von Neumann entropy scales logarithmically 
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slog" 

FIG. 6: (Color online) von Neumann entropy at half chain for 
6 — 0.3-7T and two values of D: D = — 2J (circles) and D = 
(squares). The entropy is plotted as a function of | log |r so 
that the entropy should behave as a straight line whose slope 
is c. The solid and dashed lines are the corresponding linear 
fits with slope 1.09 and 1.91 respectively. 



the disappearance of the dimer order parameter and the 
absence of magnetization. The ncmatic order parame- 
ter Q = (Sfj) — 2/3 is different from zero in this phase. 
However, this is not sufficient to signal unamibiguously 
this phase since Q is different from zero also in the dimer 
phase for D < 0. It is expected, as explained in Ref.—, 
that this phase is characterized by enhanced nematic 
correlations. This should be also true for the Critical 
A phase in the limit of large negative D. Since these 
phases are gapless, accurate DMRG calculations are com- 
putationally hard to perform. Nevertheless, we compute 
the block entropy in the Critical B phase for 9 = —0.7tt 
and D/J= —0.5 and obtain a central charge estimate: 
c = 0.99 which is compatible with the predicted result 
c = 1 valid in the limit of large negative D. This re- 
sult is in agreement with Ref in which a Bcrczinskii- 
Kosterlitz-Thouless transition is predicted (see also^-). 

IV. SUMMARY 



with I, Ref4£: 



5(0 =« log 
6 



L ttI 
-sin — ) 

7T L 



A 



(13) 



where c is the central charge of the corresponding confor- 
mal field theory and A is a non universal constant. Our 
results for S(L/2) computed at half chain are shown in 
Fig. E] for D = -2 J and for D = and 6 = 0.3tt. The 
best fits of function (fl3|) for D = —2 J and for D = give 
c = 1.09 and c = 1.91 respectively for the central charge. 

Moreover, our results are compatible with previous 
results^ which obtained at 8 = 7r/4 and D = the 
central charge c = 1.995, which is compatible with the 
expected result c = 2. On the other hand we expect, in 
the limit \D/J\ 3> 1 and D negative, the central charge to 
become c = 1, in agreement with our numerical results. 
In fact, our perturbative analysis of Sec. MI Al predicts 
the model to be described by an effective XX Z spin 1/2 
chain in the critical region. A numerical estimate for 
the central charge for this model is very close to c = 1 
as shown in Ref.—. By studying the change in the esti- 
mated central charge for fixed and as a function of D 
we find the line separating the two regions. Whether a 
crossover or a phase transition separate the two regions, 
remains an open problem. 

The boundaries of the Critical B phase are found by 



Summarizing, we have studied the complete phase di- 
agram of the spin-1 chain in the presence of a uniax- 
ial symmetry breaking field. Our study is based on the 
density matrix renomalization group (DMRG) algorithm 
where finite size effects have been carefully taked into 
account to determine the phase boundaries with accu- 
racy. Our results compare well to previous results of 
some points in the phase diagram calculated by quantum 
Monte Carlo algorithms, and agree qualitatively with the 
recent study presented ir>i£. Furthermore, our numerical 
results for large negative value of the anisotropy, coincide 
with those predicted by an effective spin- 1/2 XXZ model 
that we obtain as a limit of degenerated perturbation the- 
ory. Finally, we remark that the whole phase diagram is 
achievable with ultracold spinor gases by realizing that 
ground states of antifcrromagnctic phases can be mapped 
to the maximal energy states of ferromagnetic phases^. 
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